Quantum baker maps with controlled-NOT coupling 
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The characteristic stretching and squeezing of chaotic motion is linearized within the finite number 
of phase space domains which subdivide a classical baker map. Tensor products of such maps are 
also chaotic, but a more interesting generalized baker map arises if the stacking orders for the factor 
maps are allowed to interact. These maps are readily quantized, in such a way that the stacking 
interaction is entirely attributed to primary qubits in each map, if each j'th subsystem has Hilbert 
space dimension Dj = 2™ J . We here study the particular example of two baker maps that interact via 
a controlled-not interaction, which is a universal gate for quantum computation. Numerical evidence 
indicates that the control subspace becomes an ideal Markovian environment for the target map in 
the limit of large Hilbert space dimension. 

PACS numbers: 03.65.Yz, 05.45.Mt, 03.67.Lx 



I. INTRODUCTION 



The baker map displays the essential features of clas- 
sically chaotic motion in such a simplified form as to be 
almost a caricature. It may be described as a space-filling 
horseshoe map that linearizes the hyperbolic motion. 
Thus, essentially chaotic motion is exhibited, which can 
be followed through very simple computations. The bi- 
nary symbolic dynamics propagates vertical strips in the 
unit square onto horizontal strips. The vertical and hor- 
izontal rectangles are narrower, for longer binary codes, 
so that the primary digit specifies a half-square. Hence 
this digit is responsible for a coarse-grained description of 
the motion. The way in which horizontal strips are piled 
up by the mapping is also determined by the primary 
binary digit. 

The various quantization schemes that have been pro- 
posed for the baker map generally respect the main di- 
vision of the square into a pair of vertical rectangles, 
which split the position states into Hilbert subspaces. 
These are mapped respectively onto a pair of momentum 
subspaces according to the primary binary digit. The 
discrete and finite nature of the Hilbert space prevents 
the association of longer codes to ever thinner strips pj . 
Thus, to some extent, the primary digit is even more 
important in quantum mechanics. Initially, it may have 
seemed that the introduction of binary symbols to de- 
scribe evolving quantum systems could only serve as an 
artificial prop for the study of the semiclassical limit. 
However, the growing interest in quantum computation 
has brought the binary structure to the fore. Schack 
noted that the quantum baker could be efficiently real- 



* Electronic address: 
URL: http: //www. cbpf .br/~vallejos 



^Electronic address: psantoro@ifi.unicamp.br 
^Electronic address: < 

URL: http://www.cbpf.br/-ozorio 



vallcj os@cbpf.br 
o@cbpf.br 



ized in terms of quantum gates |2J. A three qubit Nu- 
clear Magnetic Resonance experiment was proposed la 
and then implemented (with some simplifications) [j]. 
In this context, the quantum baker map has reemerged 
as an ideal simplified model. Schack and Caves jfj have 
indeed developed an entire class of quantum baker maps, 
based on the 2 JV -dimensional Hilbert space of N qubits, 
which include the original quantizations of Balazs and 
Voros @ and Saraceno Q as special cases. While incor- 
porating the qubit structure and hence a Hilbert space 
of D = 2 N dimensions, we here adopt the Balazs- Voros- 
Saraceno map throughout, since this singles out the pri- 
mary bit most clearly. 

The object of our study is the coupling of quantum 
baker maps, initially restricted to a pair of maps. Clearly, 
the higher dimensional product-map of classical baker 
maps is also chaotic, in spite of possible added symme- 
tries. The minimal coupling that can be meaningfully 
addressed relies exclusively on the primary binary digit. 
Thus, each baker map is only sensitive to the most coarse- 
grained information available about the state of the other 
map. In other words, one can investigate the class of cou- 
pled baker maps that interact through the order in which 
the pair of horizontal rectangles are piled up. Once these 
coupled maps are quantized, this corresponds to an in- 
teraction among only the primary qubits. 

It is important at this stage to distinguish two alter- 
native forms of classical-quantum correspondence for in- 
dividual baker maps. In the original quantizations 0,0, 
the classical phase space is two-dimensional and it is the 
chaotic motion in this two-dimensional torus that corre- 
sponds to the quantum level repulsion, in keeping with 
the spectra of random matrices 0. If the quantization 
is chosen with a value of Planck's constant such that the 
Hilbert space has exactly 2 N states, then it is possible to 
reinterpret the quantum system to be the tensor product 
of N two-level systems. This would correspond to a hy- 
percube in a 2iV-dimensional phase space, but it would be 
stretching it to ascribe a classical correspondence to this 
quintessentially quantum system. We shall here keep to 
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the original quantization schemes and, hence, interpret 
each quantum baker map as corresponding to a classi- 
cal two-dimensional phase space, so that the coupling of 
two baker maps corresponds to a four-dimensional phase 
space. The four possible combinations of primary binary 
digits -(0, 0), (0, 1), (1, 0), (1, 1)- then correspond to four 
squares in the classical position space (gi, g2) and hence 
to four parallelepipeds in the unit 4-cube (gi, g2,pi,p 2 ) 
in which the classical motion is defined. 

Perhaps the simplest choice for coupling two classical 
bakers is a simple switch of the piling: Each map stacks in 
accordance to the other's primary digit. Obviously, there 
is then no change if their digits coincide, but, otherwise, 
both maps invert their piling. This is analogous to a 
mutual spin flip for the principal qubits of the quantized 
"switch baker map". Notice that this interaction does 
not break the exchange symmetry of the product map. 

As a direct model for (quantum) computing, we here 
concentrate on the "controlled-not baker map" . In this 
unsymmetric case, the control map evolves just as an or- 
dinary baker map, while it determines the stacking of the 
target map. The latter piles its rectangles as a normal 
baker when the primary digit of the control map is zero, 
but this is switched when the control digit is one. It is 
well known that the controlled-not gate is a basic element 
of quantum computing |10| . The present may be consid- 
ered as an interesting example of a system where only a 
singe element of a higher dimensional Hilbert space de- 
termines a qubit. This has been advocated as possibly 
beneficial for the stabilization of quantum computation 
against decoherence [ll|. 

It is possible to study the evolution of the entangle- 
ment of the various qubits resulting from the iteration of 
an individual quantum baker map 8] . In the case of cou- 
pled baker maps, we can interpret the primary bit of each 
component baker map as implicated in quantum compu- 
tation, while the remaining qubits in each map model 
the local environment, leading to a loss of quantum co- 
herence. Alternatively, we here consider a second baker 
map as a model for the environment coupled to the map 
that is singled out as the open quantum system. 

The presentation of the class of minimally coupled 
maps is preceeded by a brief review of the ordinary two- 
dimensional baker map and its quantization in section 
2. Section 3 then presents the coupled maps and their 
quantization with special emphasis on the controlled-not 
baker map. The evolution of bipartite entanglement of a 
product pure state is studied in section 4 by considering 
the control and the target of the controlled-not baker as 
the separate subsystems. The asymmetric nature of this 
interaction suggests that this should be an ideal Marko- 
vian system, in the appropriate limit, as far as the target 
is concerned. This is verified by the evolution of the lin- 
ear entropy. Our results are discussed in section 5. 



II. REVIEW OF CLASSICAL AND QUANTUM 
BAKER MAPS 

In this section we present the well known classical and 
quantum ingredients of baker maps. However, the alter- 
native pilings, which are usually equivalent, need to be 
made explicit when different maps are coupled. 

The classical baker transformation is an area preserv- 
ing, piecewise-linear map, b : (j>o,qo) —> (j>i,qi) of the 
unit square (periodic boundary conditions are assumed) 
defined as 

Pi = ^(Po + e ) , qi = 2g - £o ; (l) 

where eo = [2go], the integer part of 2g . This map is 
known to be uniformly hyperbolic, the stability expo- 
nent for orbits of period L being Llog2. Moreover it 
admits a useful description in terms of a complete sym- 
bolic dynamics. A one to one correspondence between 
phase space coordinates and binary sequences, 

(p,q) f-> ...e_ 2 e-i • e eie 2 ... ,£, = 0,1, (2) 

can be constructed in such a way that the action of the 
map is conjugated to a shift map. The symbols are as- 
signed as follows: is set to zero (one) when the i-th it- 
eration of (p, q) falls to the left (right) of the line q = 1/2, 
i.e. [2qi] — £j. Reciprocally, given an itinerary 

• ..e_ 2 e-i-e eie 2 ..., (3) 

the related phase point is obtained through the specially 
simple binary expansions 

oc oo 

8 = E ^St > p = E • ( 4 ) 

i=0 i=l 

Once the dynamics has been mapped to a shift on bi- 
nary sequences it is very easy to analyze the dynami- 
cal features of the map. In particular, periodic points 
are associated to infinite repetitions of finite sequences 
of symbols. 

Due to its piecewise linear nature, the baker map ad- 
mits a (mixed) generating function which is a piecewise 
bilinear form, 

W ea (pi, g ) = 2pig - e p 1 - e g , e = 0, 1 . (5) 

It is not defined on the whole space (pi,go) but on the 
classically allowed domains 

Rq = [0,1/2]® [0,1/2] and ife = [1/2, 1] ® [1/2, 1] . (6) 

Though the above generating function will be the starting 
point for quantization, it must be remembered that it 
only provides an implicit formula for each iteration of 
the classical baker map. 

With respect to the quantum map, we will follow the 
original quantization of Balazs and Voros 0, as later 
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modified by Saraceno |7j to preserve in the quantum 
map all the symmetries of its classical counterpart. In 
the mixed representation the baker's propagator can be 
written as a D x D block-matrix (D even): 

ip m \B D \q n )=^ G l' 2 G ° ) , (7) 

where position and momentum eigenvalues run on a dis- 
crete mesh with step l/D = h (h — Planck's constant), 
so that 

q n = (n+l/2)/D , p m = (m+l/2)/D , 0<n,m<D-l; 

(8) 

and Gd is the antiperiodic Fourier matrix, which trans- 
forms from the q to the p basis, 

G D = (p m \q n ) = (l/VD)e- 2mD ^ . (9) 

It will be useful to consider this matrix as the position 
representation of the Fourier operator Gd, such that 

G D \pn) = \q n ) ■ 

The mixed propagator J7J) has the standard structure 
of quantized linear symplectic maps |12| , 

(p m \B D \q n ) = l v ^Se-«hrom(p m , fc ) if e Rl 

[ otherwise . 

(10) 

In this quantization, only those transitions are allowed 
that respect the rule [2p m ] = [2q n ], a reflection of the 
classical shift property. 

To be able to iterate the quantum baker map, the 
state must be brought back to the position representa- 
tion. This is achieved by an inverse Fourier transform, 
so that the matrix (Q is multiplied by G^ 1 . 

Following Schack and Caves |5j] , we can reinterpret this 
quantum map as the evolution of N qubits if the dimen- 
sion of the Hilbert space satisfies D — 2 N . Then the 
position states can be defined as product states for the 
qubits in the basis, 

\q n ) = Id) ® lea) ® -M, (11) 
where n has the binary expansion 

N 

n = e 1 ...e N =Y,tj2 N ~ j (12) 
i=i 

and q n = (n + 1/2) /D = • ei...ejvl. The connection 
with the classical baker map is specified by the symbolic 
dynamics. Finite segments of the bi-infinite strings © 
that determine points in the unit square are made to 
correspond to orthogonal quantum states. Half of the 
position states lie in either of the two rectangles, Ro, or 
Ri defined in ©, which correspond respectively to 0, or 
1 eigenstates of the principal qubit. 

The full unitary operator for the quantum baker can 
be written explicitly as 

B D =G D [l2®G D ) 2 l (13) 



where I2 is the unit operator for the first qubit and G D , 2 
is the inverse Fourier operator on the remaining qubits. 
The operator in the square brackets preserves the first 
qubit, while evolving separately the remaining qubits, 
within each domain R Cl . It is the final Fourier operator 
that mixes the principal qubit in with the rest, because 
it acts globally on the states in both domains. This step 
is not explicit in the mixed representation J7J. 

So far we have only allowed for a single possibility in 
which to stack the rectangles in the baker transformation, 
but alternative to QJ, the classical map b' : (po,qo) — > 
{Pi,qi), 

Pi = ^(Po + 1 - eo) , qi = 2q a -e ; (14) 

has very similar properties. This variation is tantamount 
to reversing the primary classical bit and it corresponds 
to the quantum map B' , represented by the matrix: 

(p m \B D \q n }=( G ° D/2 G J/ 2 ) . (15) 

In other words, we here substitute the operator I2, which 
acted on the principal qubit, by X2 represented by the 
Pauli matrix: 

*»=(;;), as) 

so that B' D = G D [X 2 <8> G' 1 ^]. 

Of course, we are free to substitute any other unitary 
operator acting on the primary qubit, but it is only B' 
that can be interpreted classically as an equivalent alter- 
native piling of the baker map. We can still split up the 
evolution into domains equivalent to 10, though with a 
different matching of qo and p\ segments. The classi- 
cal evolution within each domain is again determined by 
classical generating functions like (J5J which become the 
exponent of the propagator. 

If we allow a finite probability for a stacking fault in 
the classical baker map, the evolution acquires a random 
component. However, if the piling order depends on the 
coarse-grained position of another baker map, the over- 
all motion will again be purely deterministic. Being that 
each baker map is thoroughly chaotic, it will be hard to 
distinguish the random motion of one of the components 
taken on its own from that of a chaotic system with an 
added stochastic variable. In the following section we 
allow the baker to choose between these alternatives, de- 
pending on its interaction with another baker. 



III. COUPLING BY THE PRIMARY QUBIT 

All possible couplings of a pair of quantum qubits are 
specified by unitary matrices acting on the basis states 
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(00, 01, 10, 11). In particular, 



h = 



10 
10 
10 
1 



(17) 



leaves both qubits invariant. This is the correct descrip- 
tion for the pair of principal qubits for two uncoupled 
quantum baker maps, B ® B. Likewise, a pair of un- 
coupled bakers, B' ® B', with the alternative stacking 
discussed in the previous section are propagated by the 
matrix 



Xo — 



( 1 

10 

10 

V i o o o 



(18) 



Perhaps an analogy may bring home the subtle nature 
of the interaction that is now introduced between the two 
subsystems. A classical analogue for these finite numbers 
of qubits could be a pair of card packs. We do not allow 
here any exchange of cards between the packs, such as 
in a chromosome crossover. Each pack is always shuffled 
separately by splitting in two and repiling, but each af- 
fects the order of the other's stacking. The interesting 
interactions are the ones which may, but do not always 
affect the piling order of the component baker map. One 
possibility is the 



SWAP 



10 
10 
10 
1 



(19) 



Here the individual bakers swap their stacking, but there 
results no change if they already shared the same digit, 
or 1. 

Perhaps, the nicest example in the context of quantum 
information is the controlled-not gate, CNOT, represented 
by the unitary matrix: 



CNOT = 



f 1 

10 

1 

\0 1 



(20) 



This is one of the universal elements for the design of 
circuits in the theory of quantum computation. Here, 
the interaction goes only one way: If the corresponding 
classical control bit were 0, the target bit would not be 
changed, propagating the full map as b, but this switches 
to b' for the control 1. In quantum mechanics the ampli- 
tudes for both choices are superposed, which entangles 
the primary control qubit to the target qubits, so that 
eventually the full set of qubits of both maps are also 
entangled. 

The general construction of interacting quantum baker 
maps is now presented through the example of the 



controlled-not interaction. Naturally, the control qubit 
now becomes the control subsystem with Hilbert space 
dimension D c , whereas the target subsystem has dimen- 
sion D t . Even if the interaction is more symmetric, D c 
will be the dimension of the Hilbert space related to the 
qubit that specifies 2x2 blocks in the above 4x4 ma- 
trices, whereas D t is related to the qubit within each 
block. The structure of the unitary operator for inter- 
acting baker maps is then 



-BCNOT = {GD c ®G Dt ) [CNOT (8) (G 



D c /2 ' 



>G- D ] /2 )}. (21) 



This is an obvious generalization of the quantum baker 
map (|13fl . The Fourier transforms become tensor prod- 
ucts of Fourier transforms and the identity for the princi- 
pal qubit is substituted by CNOT for the pair of principal 
qubits. Variations within the same structure are obtained 
by inserting other two qubit transformations instead of 
the controlled-not gate. 

The unitary matrix that describes a single step of the 
coupled maps in the mixed representation is then 



(Pm'Pm I -EcNOT I 



/M0 

M 

M 

V M 



with 



M = G 



DJ2 



G 



D t /2 ■ 



(22) 



(23) 



Each of the four nonzero blocks in the matrix above is 
a D c /2 x D t /2 (inverse) Fourier transform for the states 
corresponding to four classical domains, namely products 
of the rectangles similar to those defined by @. Within 
each of these 4-dimensional parallelepipeds, the classical 
evolution is just the linear hyperbolic motion specified 
by generating functions that are the sum of those for 
the baker maps b, i.e. @, or b'. Just as with a single 
baker map, it is the final D c x D t Fourier transformation, 
returning to the position representation, that effects an 
interaction of the principal qubits with all the others. 
But it should be noticed that there is no interaction be- 
tween the secondary qubits of the control Hilbert space 
and those of the target prior to the final Fourier trans- 
formation. 

The classical correspondence for the interacting quan- 
tum baker maps only makes sense if both the dimensions 
D c and D t are large. However, it is interesting to con- 
sider the limit where the control subsystem is a single 
qubit, i.e. D c = 2. In that case, the full map is just 



-BCNOT = {G Da ® Gp t )[CNOT (g> G 



D t /2\ 



(24) 



where Gd c is a single qubit gate (had we used the peri- 
odic Fourier transform, it would be the Hadamard gate) . 
In a way, we have here removed the internal degree of 
freedom of the control, which now only functions through 
its effect on the target. This is somewhat similar to the 
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FIG. 1: Distribution of level spacings -normalized to unit 
average- for the CNOT baker map (big dots). The dimensions 
are 32 and 16 for the control and target bakers, respectively. 
We also display the predictions of random matrix theory (dot- 
ted line: GOE, dashed: GUE) and the Poisson distribution 
(full line). 



generalizations of the quantum baker map proposed by 
Schack and Caves |5| , which single out more than a single 
qubit. In the case of a pair of qubits, their map becomes 

B N , 2 = i 2 ®[G~^ t {\ 2 ®G Dt )]S, (25) 

where we have identified D t = 2 N ~ 2 , while S permutes 
the pair of principal qubits. The resemblance becomes 
stronger if we substitute the CNOT operator in the inter- 
acting bakers by 12<8>12' This comparison shows that our 
general construction of interacting maps allows for the in- 
vestigation of a much richer range of dynamics than that 
displayed by previous generalizations of quantum baker 
maps, even when the internal structure is removed from 
one of the maps. 

In spite of the fact that the quantum controlled-not 
baker corresponds to a classically chaotic map, its spec- 
tral statistics does not comply with one of the standard 
ensembles of random matrix theory 0, as we see in 
Fig. n This discrepancy is shared by many previous ex- 
amples 0, 0. In most cases, hidden symmetries, or 
arithmetic anomalies have been identified, which distin- 
guish the system from the generic ensemble |l4j . 

The classical controlled-not baker bcNOT : (Po,<7o) —> 
(.Pij<7i)j corresponding to the unitary operator (|21|l . is 

Pi = \(Po + t c ), (26) 

<ll = 2 9 g - e c , (27) 

P\ = ^(pS + 4 + eo-2eSe§), (28) 

q\ = 2ql - el . (29) 



Thus the classical orbits depend on the binary digits, e§ 
and 6q in an unsymmetric way. The orbit of the control 
subsystem is independent of the behaviour of the target, 
while the evolution of the principal binary digit of the 
control, that results from its internal chaotic motion, af- 
fects the target. This influence is as random as a perfect 
coin toss, so the control can be considered as a perfect 
Markovian environment for the target. 

It is important to distinguish the classical trajectories 
of the full 4-dimcnsional map, from its projections into 
the 2-dimensional phase space rectangles that describe 
each of the subsystems. The full trajectory is determin- 
istic, which is also the case for its projection onto the 
control phase space. It is the projection onto the tar- 
get phase space that acquires an important random com- 
ponent, when knowledge of the control orbit is deleted. 
Viewed from the control phase space, we find an infinite 
number of trajectories of the full system that project onto 
the same control orbit. 

Viewed within quantum mechanics, the relation be- 
tween the subsystems becomes less unsymmetrical. Their 
entanglement is usually measured with reference to the 
reduced density matrix of either subsystem, that is p Cl or 
Pt, such as the purity, tr p 2 , or, equivalently, the linear 
entropy [Icj . 

S L = 1 - tr p 2 . (30) 

The effect of the entanglement is an increase of linear 
entropy from that of an initial product state, for which 
Sl = 0. These measures are necessarily the same for 
both subsystems, no matter how lopsided the interaction 
[To|. Nonetheless, it is legitimate to consider the target 
subsystem of the controlled-not baker interaction as an 
ideal candidate for a quantum Markovian system. It is 
this picture of the control as forming an environment 
for the target that shall be studied below, rather than 
the alternative investigation of the entanglement of the 
principal qubits, with local environments taken as the 
set of remaining qubits of the controlled-not baker. The 
latter is the line taken in 

IV. EVOLUTION OF ENTANGLEMENT: THE 
MARKOVIAN LIMIT 

We considered the unitary evolution of the quantum 
controlled-not baker map, given initial product states, 
\ij}) c \ip) t . Both of these factor states were chosen ran- 
domly with respect to the unitarily invariant Haar mea- 
sure for pure states within both the control subsystem 
and the target subsystem . Figure |2] show the evo- 
lution of the linear reduced entropy for either subsystem 
as a function of the discrete time, i.e. the number of it- 
erations. In each of these figures the dimension of the 
target space was kept fixed at D t = 16, that is, four 
qubits. However, the size of the environmental control 
subspace grows from D c = 4 in (a) to D c = 256 in (d), 
i.e. eight qubits. In all cases the linear entropy rises 
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FIG. 2: Linear entropy as a function of the number of itera- 
tions of the CNOT baker. The dimension of the target space 
is fixed at D t = 16, but the size of the control subspace is 
increasing: D c = 4 (a), 16 (b), 64 (c), 256 (d). In each case 
we considered 50 initial product states chosen according to 
the (product) Haar measure. Also shown are the entropies of 
corresponding random states (horizontal lines). 



FIG. 3: Markovian limit for the target map of the CNOT baker. 
Shown is the linear entropy as a function of the number of 
iterations of the superoperator of Eq. 132H . The dimension 
of the target space is Dt = 16. We chose 50 initial product 
states according to the Haar measure. 



sharply during the first five iterations, or so, and then 
oscillates around a stable average value, which depends 
on the choice of initial state. The asymptotic values for 
the entropy can be compared with those of random states 
. The average value of the linear reduced entropy for 
the Haar ensemble appropriate to the full Hilbert space 
is less than maximal 



ppropi 

urn 



(Sr.) = 1 - 



D c + Dt 
D c D t + 1 



(31) 



In each of the above figures the Haar average is seen to lie 
within the fluctuations of the evolution of the individual 
states. This means that the entangling power |l5j of the 
CNOT baker is very close to that of a random operator 
|17|. Another striking effect of increasing the dimension 
D c is the decrease in the amplitude of the fluctuations 
in linear entropy about their mean. This is in line with 
the decrease of the variance of the purity in the Haar en- 
semble as the dimension of the subsystems are increased 

Let us now compare this behaviour of the linear en- 
tropy of the target subsystem with that of the Markovian 
system: 



(32) 



That is, we average over the two possible stackings of the 
quantum baker map within the standard formalism of 
Kraus superoperators ^8|. This Markovian evolution is 
shown in Figure[3]for a set of random initial target states 
as in Fig. [21 so that for D t = 16. Evidently, the closest 
correspondence is verified with Fig. |2Id) , which has the 



largest control space. Even though the target only inter- 
acts directly with the principal qubit of the control space, 
the internal motion within the latter is required to wash 
away information concerning the evolution of the target. 

Since the Markovian approximation presupposes com- 
plete randomness for the initial state of the environment, 
it is more appropriate to compare this irreversible evo- 
lution with the unitary evolution of a density operator 
defined as a product of a pure target state with a mixed 
control state. The latter is chosen to be proportional to 
the identity operator, i.e. it has equal weight for all basis 
state projectors in the control subsystem. Figure 0] dis- 
plays the growth of the linear reduced entropy for a sin- 
gle initial state chosen randomly for various dimensions 
of the control subspace, D c . This is strong evidence that 
the unitary evolution of the linear entropy converges onto 
the Markovian evolution in the limit where D c — > oo. 

Most of the eigenvalues of the Markovian superopera- 
tor that acts on pt (0) in (|32J) lie within the unit circle, as 
shown in Figure|S] In time, only the projection of the ini- 
tial density operator with A = 1 remains. However, this 
eigenvalue is doubly degenerate, corresponding to both 
the identity operator and the reflection operator, defined 
by R\q m ) = |1 — q m ) (the reflection is a symmetry of the 
individual baker maps Q). The fact that different initial 
density operators also differ in their projections onto this 
pair of states accounts for the spread in the asymptotic 
values of the respective linear entropies. 
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FIG. 4: Comparison of Markovian vs. unitary evolution, (i) 
A single randomly-chosen initial state of the target system 
was evolved with the Markovian Eq. 1321 . (ii) The same tar- 
get state, tensored with the identity of the environment, was 
evolved with the full unitary dynamics. In both cases we plot 
the linear entropies as a function of time. The Markovian case 
corresponds to the dotted line. Full lines represent the uni- 
tary evolution for environment sizes D t = 8, 16, 32, 64. The 
dimension of the target space is fixed at Dt = 16. 
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FIG. 5: Eigenvalues of the Markov superoperator. The di- 
mension of the target is D t = 16. 



V. DISCUSSION 

The general idea of a baker map is that "vertical" 
rectangular partitions of the unit square are stretched 
horizontally and squeezed vertically and then stacked. 
The generalization to partitions of 2L-dimensional hy- 



percubes is obvious, so that these then model a system 
with L degrees of freedom. We have here shown that 
the different stacking orders allowed by higher dimen- 
sional baker maps may be interpreted as resulting from 
the various interactions of the principal qubits for each 
subsystem, if each of these is quantized to have a finite 
Hilbert space of dimension Dj = 2 nj . As a first explo- 
ration of such systems we have chosen the controlled-not 
baker map, where the interaction of the principal qubits 
is the universal CNOT gate for quantum computation. We 
have discussed other possibilities for L = 2 and there are 
many more for increasing L. 

It is also possible to couple a baker map to some other 
simple unitary evolution. In the case of Ermann, Paz and 
Saraceno , the latter corresponds to a pair of transla- 
tions. Then the control map decides the direction of the 
displacement, rather than a stacking order for the target 
map. Such a quantum random walk does not belong to 
the class of generalized baker maps, because the target 
is not a chaotic system. 

A further more elaborate possibility for interaction in- 
volves a rotation of the four parallelepipeds into which 
the primary digits of both maps divide a four dimen- 
sional phase space. Choosing this to be of 7r/4, shuffles 
the parallelepipeds, so that the doubled binar y d escrip- 
tion is respected. This system is studied in 20]. The 
interesting point is that the equilibrium and all the clas- 
sical periodic orbits become "loxodromic" , i.e. the posi- 
tions spiral outwards, while the momenta spiral inwards, 
because the eigenvalues of the stability matrix are gen- 
eral complex numbers. So far, not much effort has been 
made towards quantizing loxodromic motion, but |21| is 
an exception. 

Unlike the case of the various quantizations of a single 
baker proposed by Scott and Caves jl] , there is a natural 
partition into the separate Hilbert spaces for interacting 
baker maps, with perfect classical correspondence. We 
can analyze the entropy evolution for the reduced density 
matrix of a single subsystem for each form of coupling. 
The other clear alternative is to trace out the density 
matrix of the secondary qubits for each component baker, 
so as to obtain a mixed state for each principal qubit. 
Then their entanglement is uniquely determined by the 
concurrence p^ . 

The remarkable facility for computing numerically the 
evolution of both quantum and classical baker maps ex- 
tends to their higher dimensional couplings. Here we 
have provided numerical evidence that the generalized 
controlled-not baker map may be an ideal model for 
Markovian evolution of the target map, once the envi- 
ronmental control map is traced away. It is interesting 
to note that the interaction with the environment is not 
weak in any usual sense. 
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